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In this work we focus on the study of phase separation in the zero-bandwidth extended Hub¬ 
bard with nearest-neighbors intersite Ising-like magnetic interactions J and on-site Coulomb in¬ 
teractions U. The system has been analyzed by means of Monte Carlo simulations (in the grand 
canonical ensemble) on two dimensional square lattice (with N = L x L = 400 sites) and the re¬ 
sults for U /(4J) = 2 as a function of chemical potential and electron concentration have been 
obtained. Depending on the values of interaction parameters the system exhibits homogeneous 
(anti-)ferromagnetic (AF) or non-ordered (NO) phase as well as phase separation PS:AF/NO state. 
Transitions between homogeneous phases (i.e. AF-NO transitions) can be of first or second order 
and the tricritical point is also present on the phase diagrams. The electron compressibility K is an 
indicator of the phase separation and that quantity is of particular interest of this paper. 

PACS numbers: 

71.10. Fd — Lattice fermion models (Hubbard model, etc.) 

71.50.-b — General theory and models of magnetic ordering 
75.30.Fv — Spin-density waves 

64.75.Gh — Phase separation and segregation in model systems (hard spheres, Lennard-Jones, etc.), 

71.10. Hf— Non-Fermi-liquid ground states, electron phase diagrams and phase transitions in model systems 
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I. INTRODUCTION 

Magnetic insulators are a class of materials realized in 
many various compounds. An example of them are ma¬ 
terials known as transition metal cluster compounds with 
general formula AM, jA'g, where A —trivalence metal, 
M —transition metal, X —chalcogenide m- Phase sep¬ 
arations can occur in magnetic insulators in various cir¬ 
cumstances and their theoretical understanding is very 
current topic. Moreover, instabilities such as stripe for¬ 
mation (as well as charge order) can occur in high-T c su¬ 
perconductors (e.g. in cuprates) [4j. A simplified model 
to describe behavior of such materials is the extended 
Hubbard model with intersite magnetic interactions 13 
udi- In this work we study the zero-bandwidth limit of 
the extended Hubbard model. The Hamiltonian of this 
model has the form: 

H = U nifhn + 2 J ^ s* - n ^ fij (1) 
* (*b> » 

where U is the on-site density-density interaction, J is 
2 -component of the intersite magnetic exchange interac¬ 
tion, and fi is chemical potential. The interactions are ef¬ 
fective model parameters and are assumed to include all 
the possible contributions and renormalizations. 
restricts the summation to nearest neighbors (indepen¬ 
dently) . hi = hif + hii is total electron number on site 
i and sf = (l/2)(hn — fin,) is 2 -component of total spin 
at i site. hi a = cf a Ci a is electron number with spin a on 
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site i, where cf a and Ci a denote the creation and annihi¬ 
lation operators, respectively, of an electron with spin <j 
(a =t, i) at the site i. The electron concentration n is de¬ 
fined as n = (1/iV) JA(fij), where N is the total number 
of sites. 

The model studied exhibits two symmetries: (i) the 
symmetry between J > 0 (antiferromagnetic) and J < 0 
(ferromagnetic) cases and (ii) the electron-hole sym¬ 
metry. Because of these symmetries only analyses for 
0 < n < 1 and J > 0 have been performed. 

We have used the Monte Carlo (MC) simulations to 
analyze the system. Simulations have been done using 
Hamiltonian described above on two dimensional square 
lattice with N = L x L sites in the grand canonical en¬ 
semble, which allows us to obtain e.g. chemical poten¬ 
tial dependence of electron concentration curves — n(/z). 
The Monte Carlo algorithm used in this analysis con¬ 
sists of three steps: (i) creation, (ii) destruction, and (iii) 
moving of particle, all of them with appropriate proba¬ 
bility P ~ exp (/S.E/(ksT)) |lltiT4| . It is worth noting, 
that for constant values of concentration a simpler algo¬ 
rithm with only step (iii) — „move” would be sufficient. 
However, addition of the grand canonical parts (creation 
and destruction) allows one for more detailed analysis in 
full range of chemical potential and concentration. Un¬ 
fortunately, addition of chemical potential term in the 
Hamiltonian prevents us from implementing cluster up¬ 
dates algorithm, so only local updates m are used here. 
The details of the algorithm used can be found in EEHU]. 

The exact ground state (T = 0) results for this model 
have been found in the case of a d = 1 chain OH HU using 
the Green function formalism as well as for 2 < d < +oo 
case mm • The rigorous results for finite tempera- 
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FIG. 1. The phase diagrams of the model for C//4J = 2: (a) fcsT/4J vs. /i/4 J and (b) fcflT/4J vs. n (L = 20). T denotes a 
tricritical point. On panel (a) “heating” and “cooling” label the boundaries obtained by simulation performed with increasing 
and decreasing temperature, whereas “average” is the average of these two results. They differ if the AF-NO transition is 
first-order. 


tures T > 0 have been also obtained |T7| [20] for d = 1 
chain (an absence of long-range order at T > 0). Within 
the variational approach (with on-site U term treated ex¬ 
actly and mean field decoupling of intersite term J) the 
model has been analyzed for half-filing (n = 1) [ST) 122 
as well as for arbitrary electron concentration 0 < n < 2 
[25] (these results are rigorous in the limit of infinite di¬ 
mensions d —> Too). Our preliminary Monte Carlo (MC) 
results have been presented in puna for L = 10 and on¬ 
site repulsion: [//(4J) = 1,10 (corresponding to rather 
weak and strong coupling, respectively) [13] as well as 
for L = 20 and [7/(4 J) = 1 |I3] . 

In the present paper we investigate in details the phase 
diagram and thermodynamic properties of the model 
for arbitrary electron concentration n < 1 and arbitrary 
chemical potential p < 0 (ji = fi — U /2) in the whole 
range of temperatures for 17/ (4J) = 2 and L = 20. In 
particular, we focus on a behavior of an electron com¬ 
pressibility. The corresponding results for n > 1 (p > 0) 
are obvious because of the electron-hole symmetry of the 
model on alternate lattices mentioned previously. 


termined using hysteresis analysis 031 • In Fi g- 0 a ) 
“heating” and “cooling” label the boundaries obtained 
by simulation performed with increasing and decreasing 
temperature whereas “average” is the average of these 
two results. They differ if the AF-NO transition is 
first-order. Details of this method can be found in nn. 

With simulations done for fixed /t and 
ksT/AJ vs. p/AJ, it is possible to obtain phase 
diagrams as a function of n (shown in Fig. 0 b )) b y 
determining electron density above (n_) and below 
(n_(_) the AF-NO phase transition (for fixed p). The 
first-order AF-NO boundary for fixed p splits into 
two boundaries (i.e. PS-AF and PS-NO) for fixed n. 
At sufficiently low temperatures, i.e. below T-point, 
a phase separated (PS: AF/NO) state occurs. The PS 
state is a coexistence of two (AF and NO) homogeneous 
phases. At higher temperatures (i.e. above T-point) the 
AF-NO transition is second-order one. 

An objective indicator of a PS state existence is the 
evolution of the compressibility K of the system mm- 
For a system with variable number of particles it can be 
defined as 


II. RESULTS AND DISCUSSION (1//(4J) = 2) 

Finite temperature phase diagrams for this model were 
obtained using MC simulations for U / (4J) = 2 (and 
L = 20) as a function of p/AJ and n are presented in 
Fig. 0a) and Fig. 0b), respectively. 

The behavior of this system for fixed p is rather 
simple (Fig. 0a)), with both first-order (below 
T-point) and second-order (above T-point) phase 
transitions separating non-ordered (NO) and antiferro¬ 
magnetic (AF) phases with tricirtical point T located 
at hgT/AJ = 0.205 ± 0.003, p/4J = -1.405 ± 0.007 
(n ~ 0.58). The location of T-point has been de¬ 


— =n 2 ( — ^ 

K \dn J T ,u,J 


( 2 ) 


From this definition it follows that at a fixed /t the num¬ 
ber of particles in an open system can fluctuate freely 
(precisely, in some define range) when K —» Too. Such 
a behavior is connected with an occurrence of the PS 
states in define range of n. At the same constant total 
free energy of the system the number of domains as well 
as their distribution can change. Hence, the phase sepa¬ 
ration states are „highly unstable” in that sense that they 
are subjected to continuous fluctuations of local density 
(but the total density n is constant). 
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FIG. 2. Compressibility K as a function of temperature T for 
constant values of electron concentration n (as labeled) for 
L = 20. For n > 0.58 there are fluctuations of K associated 
with second-order AF-NO phase transitions. 


concentrations than those corresponding to the T-point 
(n > 0.58) there are compressibility fluctuations related 
to second order AF-NO phase transition as shown on 
Fig- |5| at higher temperatures. They are significantly 
smaller than the fluctuations close to the boundaries be¬ 
low which the PS:AF/NO state occurs. 

Fig. [3] presents a map of compressibility on fcgT/4 J-n 
plane together with the phase boundaries derived in 
Fig . [ljb) plotted for the comparison. An increase of com¬ 
pressibility close to the second-order AF-NO boundary is 
clearly seen. The boundary between homogeneous phases 
and phase separation state is also visible, with the com¬ 
pressibility close to the boundary being at least an order 
of magnitude greater than those inside the homogeneous 
phases. As it was said earlier, in the case of phase sepa¬ 
ration occurrence K —> 00 , so no points are shown inside 
the region of the PS state occurrence in Fig. [3] 

III. FINAL COMMENTS 
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Notice that the results presented are in good qualita¬ 
tive agreement with mean field calculations using vari¬ 
ational approach presented in HU |2IH2Sj. When com¬ 
paring these results one should keep in mind differences 
between these two methods, as the VA is exact only for 
infinite dimensions d — > 00 . The drawback of Monte 
Carlo simulations is long thermalization time, which pre¬ 
vents us from obtaining results for the ground state and 
very low temperatures, as in these conditions electrons 
have very small probability of escaping local energy min¬ 
ima. Behavior of the model considered in the case of 
finite band (t ^ 0) is very interesting and mostly open 
problem in the general case mm- 
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FIG. 3. A map of compressibility Kn 2 on ksT/AJ-n plane 
(for L = 20). Solid black curves indicate the phase boundaries 
derived in Fig. 0b) plotted for the comparison. 

In Fig. [2]the compressibility Kn 2 is plotted versus tem¬ 
perature hsT/AJ for several constant values of n. As is 
clearly seen, close to the boundaries of the PS:AF/NO 
state occurrence plotted in Fig [TJb), the value of com¬ 
pressibility abruptly increase and K —► +00 at transi¬ 
tion temperature, indicating an existence of the phase 
separation state. At higher temperatures and for larger 
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